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Recent lattice results on QCD at finite temperatures and densities are reviewed. Two 
CD ' new and independent techniques give compatible results for physical quantities. The 
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^ ' phase line separating the hadronic and quark-gluon plasma phases, the critical endpoint 
and the equation of state are discussed. 

^ , 1. Introduction 
> 

^ I QCD at finite temperatures (T) and/or chemical potentials (/i) is of fundamental im- 

portance, since it describes particle physics in the early universe, in neutron stars and in 
heavy ion collisions. According to the standard picture, at high T and/or high density 
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CN ■ there is a transition from a state dominated by hadrons to a state dominated by partons. 
^ . The expression "transition" (which occurs at the transition temperature Tc) is used for 
i first/second order phase transitions and crossovers. (Observables change rapidly during 
a crossover, but no singularities appear.) Extensive experimental work has been done 
with heavy ion collisions at CERN and Brookhaven to explore the /x-T phase boundary 
^ ■ at relatively small values. At large a rich phase structure is conjectured [ |], 0, ^ If. 

There are well established nonperturbative lattice techniques to study this transition 
^ ! at vanishing density. At fi=0 we have fairly good description of the transition (e.g. the 
^ [ order of the transition as a function of the quark masses or the the equation of state as a 
function of T). For recent reviews see the summaries of the lattice conferences [ Q or the 
review talk on finite T lattice QCD at this conference [ ^. 

Our knowledge is far more limited at nonvanishing /i. Due to the sign problem (oscil- 
lating signs lead to cancellation in results, a phenomenon which appears in many fields of 
physics) nothing could have been said for almost 20 years about the experimentally impor- 
tant case at nonvanishing densities. In the last year new, and for the first time successful 
approaches appeared and physically relevant results were obtained [ 0, §|, ^, |T^, 0, |T2], |T3| . 

The aim of this review is to give a self-contained picture on the lattice approach at 
fi^O for real QCD (QCD-like models, such as SU(2), random matrix or NJL models are 
not discussed). In Section 2 the qualitative features of the phase diagram are summarised 
both at /i=0 and fJ^^O. Section 3 briefly presents the lattice formulation and shows the 
associated sign problem. The main emphasis is put to the origin of the problem and 
technical details are not discussed. In Section 4 the two new techniques (the overlap 
improving multi-parameter reweighting [ |^ and the analytic continuation [[T^) are pre- 
sented. Readers who are not interested in the origin of the sign problem and in the new 
techniques should skip Section 3 and 4. In Section 5 results of the Budapest group [ 
^ H, |13[ are listed (the phase line, the location of the critical endpoint and the equation 
of state at /^t^O). They are obtained by the direct application of the overlap improving 
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Figure 1. QCD phase diagrams. Left panel: The phase diagram on the mud versus 
plane at fi=0. The left lower and the right upper corners are first order phase transition 
regions. Thick solid lines (A-E,E-F and H-G) indicate second order phase transitions. 
The physical point is most probably in the crossover region. Right panel: The phase 
diagram on the T versus /i plane. For small T and fi the system is in the hadronic 
phase, for large T it is a quark- gluon plasma, whereas for small T and large fi colour 
superconductivity is conjectured. The lower curve shows the nf=2 case. In the small 
fi region the dashed line indicates a second order phase transition which is connected 
by a tricritical point P to a first order phase transition line. The upper curve shows 
the n/=2+l case, with physical quark masses. The second order phase transition line 
disappears, instead of it we are faced with a crossover, illustrated by the dotted region. 
This region is connected to a first order phase transition line by a critical endpoint E. 



multi-parameter reweighting. Section 6 gives the findings [|^, |r2| of the Bielefeld- Swansea 
group (the phase line, the equation of state at /i~0 and T^Tc and the response of the 
critical endpoint to the change of fi). They also use the multi-parameter reweighting 
technique; however, the /i-dependence of the determinant is approximated by a Taylor 
expansion. Section 7 discusses the phase line obtained by analytic continuation [ |10|, [Til] . 
Section 8 concludes. 



2. Qualitative features of the phase diagram 

Our knowledge on the phase diagram of QCD at fi=0 and fiy^O is summarized on 
Figure |l|. Some ingredients are rigorous lattice results, others are indications from models. 

What are the characteristics of the fi-T phase diagrams, relevant for heavy ion colli- 
sions? (See the right panel of Figure |1|.) One of the most interesting features of the phase 
diagram is a critical endpoint E connecting the first order phase transition line with the 
crossover region, which separates the low T hadronic and high T quark-gluon plasma 
phases. It is a long-standing open question, whether such a critical point exists on the 
fi-T plane, and particularly how to predict theoretically its location [ 

Let us discuss first the fi=0 case (see the left panel of Figure |l|; note, that now we study 
only those features of the figure, which are relevant for the critical endpoint at fi^O, for 
more details see [ ^ and references therein). Universal arguments [ |15[ and lattice results 



|] indicate that in a hypothetical QCD with a strange (s) quark mass {rus) as small as 
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the up (u) and down (d) quark masses (mu^d) there would be a first order phase transition 
at finite T (point in the E-B-F region). The nf=2 case with small (but nonvanishing) 
u/d quark masses but with = oo there would be no phase transition only an analytical 
crossover. This means that between the two extremes there is a critical strange mass (m^) 
at which one has a second order finite T phase transition. nj=2+l lattice results with 
two light quarks and rris around the Tc indicated that is about half of the physical m<j. 
Thus, in the real world we probably have a crossover. (Clearly, more work is needed to 
approach the chiral and continuum limits [ ^.) 

At nonvanishing n, arguments based on a variety of models predict a first order finite T 
phase transition at large /i. Combining the /i = and large /i informations an interesting 
picture emerges on the /i-T plane. For the physical rus the first order phase transitions at 
large /i should be connected with the crossover on the /x = axis. This suggests that the 
phase diagram features a critical endpoint E (with chemical potential fiE and temperature 
Te), at which the line of first order phase transitions {fi > fiE and T < Te) ends [ |14|. At 



this point the phase transition is of second order and long wavelength fiuctuations appear, 
which results in characteristic experimental consequences, similar to critical opalescence. 
Passing close enough to {^e^Te) one expects simultaneous appearance of signatures, which 
exhibit nonmonotonic dependence on the control parameters [ |1^], since one can miss the 
critical point on either of two sides. The location of this endpoint is an unambiguous, 
nonperturbative prediction of the QCD Lagrangian. Unfortunately, no ah initio, lattice 



work was done earlier to locate the endpoint. Only results from models [[l^ were available. 

The goal of present /it^O lattice studies is to determine the phase diagram, to locate 
the endpoint and also to calculate the equation of state. 



3. Lattice QCD and the sign problem at /xt^O 
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Figure 2. Illustration of Euclidean lattice QCD. Anticommuting quark fields are defined 
on the lattice sites. Gauge fields are SU(3) matrices and used as link variables: Uy{n). 
The product of four links around a plaquette leads to the standard gauge action. 



The continuum QCD Lagrangian in Euclidean spacetime is £ = — ;|-F^^F'''^'^ + ^/'(D^7^ + 
m)il). The first term gives the gauge, the second one the fermionic contribution. In the 
discretised lattice formulation the anti-commuting %l){n) quark fields live on the sites {n) of 
the lattice, whereas the gluon fields are used as link {U^{n)) and as plaquette {P^y{n)) 
variables 



U^{n) = exp {igs 



dx'^Al{x')Xa/2), 



U^{n)U,{n + eM{n + eMin) (1) 
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Here n represents the direction and the corresponding unit vector. Similarly to the 
continuum formulation the action S = Sg + Sf consists of the pure gluonic and the 
fermionic parts. The gluonic part is written with the help of the plaquettes: Sg = 
Q/ds ' ^n,fi,u [1 ~ R-e(TrP^j,(n))], where gs is the gauge coupling. The fermionic part on 
the lattice needs a differencing scheme for quarks. In the naive, noninteractive case one 
obtains: il){x)'^^d^il){x) — ^ ilj{n)'~f^{tp{n + e^) — ipln — e^)). The interacting case has 
instead of d^, thus the gauge field is also included: ip{x)'y^D^ilj{x) — >• ip{n)'y'^U^{n)ip{n + 
e^j) + ... (fermion doublers do not play any role in the sign problem, therefore we do 
not discuss them). Fermion fields appear only in bilinear expressions, thus we can write 
Sf = ip{n)Mnm4'{''T^)j where M is the [/-dependent fermion matrix. Note, that the number 
of raws or coloumns of M is proportional to the lattice volume. In our case the matrix M 
is sparse, only diagonal (oc m) and next-neighbour (oc U^) elements are nonvanishing. 

Our system can be described by the Euclidean partition function. It is given by inte- 
grating the Boltzmann weights over the gauge and fermionic fields. The action is bilinear 
in the fermionic variables, thus this part of the integral can be calculated explicitely. 

Z = y VUVipV'ilje-'^'-'^f= J VUe-^' det M{U). (2) 

The remaining integral over the U fields is calculated by stochastic methods. A canonical 
ensemble of field configurations is generated by Monte-Carlo algorithms. Observables are 
obtained as averages over the field configurations. The intrinsic feature of these techniques 
is importance sampling, thus we sample only the most important configurations and these 
individual configurations have equal weights in the averages. 

We use the most straightforward technique, the Metropolis algorithm as an illustration 
(the Metropolis method is very CPU-intensive, thus usually faster but more complicated 
techniques are used). The basic step of the method is a stochastic, reversible modification 
of a link variable: U^{n) U'^{n). It is accepted or rejected. The Metropolis condition 

P{U U') = min [1, exp{Sg{U) - Sg{U')) det M{U')/ det M{U)] (3) 

gives the probability to accept the new U' field variable. It is easy to show, that by 
sweeping through the whole lattice many times, we reach the canonical distribution, and 
physical observables can be determined. Note, that for real gauge fields det M{U) is real, 
thus eq. (0) really has a probability interpretation. This feature is essential not only for 
the Metropolis algorithm, but for any importance sampling based method. 

Our primary goal is to understand QCD at finite densities. As usual, finite densities 
can be studied by introducing a chemical potential. One adds the following term to the 
continuum action: /i^/'(x)74?/'(x). In this formula n acts as a fourth component of an 
imaginary, constant vector potential [|T^, |T8|. As we have seen real gauge fields result in 



real detM(?7); however, the inclusion of // (thus a constant imaginary gauge field) gives 
complex det M{U). This spoils the probability interpretation of eq. (^ and any other 
importance sampling based method. Instead of an ensemble of equally important config- 
urations (importance sampling) we might have configurations with complex Boltzmann 
weights. These complex weights with oscillating real parts largely cancel each other in 
observables. In the literature this phenomenon is referred to as the "sign problem" . 
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Figure 3. Illustration of the multi-parameter reweighting. The horizontal arrow shows 
the Glasgow, the curved one the new, multi-parameter reweighting. The dashed transition 
line separates the different phases. The Glasgow method starts with an ensemble in the 
low T, hadronic phase and uses the weighting factors of the configurations at the same 
T but at fi^O. Clearly, the hadronic configurations do dot overlap with the transition 
ones. The new method starts e.g. at the transition point at fi=0 and changes T and fi 
simultaneously, keeping the system on the transition line. The overlap is much better. 



4. New lattice techniques at finite chemical potential 

The overlap improving multi-parameter reweighting [ |^ opened the possibility to study 
lattice QCD at nonzero T and /x. First one produces an ensemble of QCD configurations 
at /i=0 and at Tt^O. Then the Ferrenberg-Swendsen type reweighting factors [ ^ of 
these configurations are determined at /it^ and at a lowered T. The idea can be easily 
expressed in terms of the partition function 

(3)= J VUe-^o^^^^^ det M(/i, m, U) = (4) 



det M{fi,m, U) 
detM(/i = 0,m, U) 



where Sg{U) is the action of the gluonic field, [3 = fixes the coupling of the strong 
interactions (g). Note that for a given lattice T is an increasing function of /?. The quark 
mass parameter is m and det M comes from the integration over the quark fields (see 
eq.(^). At nonzero fi one gets a complex detM which has no probability interpretation, 
thus it spoils any importance sampling. Thus, the first line of eq. at n^O, is rewritten 
in a way that the first part of the second line is used as an integration measure (at /i=0, 
for which importance sampling works) and the remaining part in the curly bracket is 
measured on each configuration and interpreted as a weight factor {w{j3, fi,m,U)}. (For 
rifj^A fractional powers of detM is needed. This complication can be solved [ |[.) 

The reweighting is performed along the best weight lines on the plane (or equiva- 
lently on the /z-T plane). One such line is the transition line. The best weight lines are 
determined by minimising the spread of logw. The technique works for T at, below and 
above T^. Using the above weights any observable O can be determined at /xt^O 



(5) 
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The method is illustrated on Fig. |^. Using multi-parameter reweighting one simulta- 
neously changes T and fi. A much better overlap can be obtained by the multi-parameter 
reweighting than by the single /i-reweighting Glasgow- method The Glasgow- method 
reweights pure hadronic configurations to transition ones. By the new technique transi- 
tion (or hadronic/QGP) configurations are reweighted to transition (or hadronic/QGP) 
ones. Since the original ensemble is collected at /i=0 one does not expect that even the 
new technique is able to describe the physics of the large /i region with e.g. colour su- 
perconductivity. Fortunately, the typical values at present heavy ion accelerators are 
smaller than the covered region. 

An alternative approach [ ^ uses Taylor expansion of eqs. (H,|D as a function of /i (or 
m) and hence estimates the derivatives of various quantities with respect to fi (or m). 
Thus, for the case of a Taylor expansion in n one has 

i«M(o) j " S ^ — ojr^ — = S • 

Instead of using the explicit form of the determinants in eq. (^) the Bielefeld-Swansea 
group used the derivatives in yu. Compared to the explicit calculation of the determinants, 
the Taylor technique needs less CPU-time (derivatives can be estimated stochastically); 
however, only valid for somewhat smaller fi values than the full technique. 

The other promising and absolutely independent approach is the analytic continuation 
from imaginary to real chemical potentials [ |10|. One computes first the critical line 
for imaginary chemical potential. For these /x values there is no sign problem, therefore 
Monte-Carlo simulations based on importance sampling can be carried out. One can check 
the convergence of the Taylor expansion of the critical line, thus f3c as a function of Im(/i). 
The convergence seems to be surprisingly fast for the whole range of chemical potentials 
accessible to this method. Analytic continuation of the Taylor series then reduces to 
simply fiipping the sign of the appropriate terms. Finally, the infinite volume limit has 
to be taken from the continued results at real /i. Using this technique the phase line 
separating the hadronic and quark-gluon plasma phases were determined [ [111. The 
order of the transition can then be determined in a (V, /i)-range where the truncation 
error of the series is smaller than finite size scaling effects. 



5. Overlap improving multi-parameter reweighting: direct approach 

This section summarizes some results of the Budapest group [ ^ |^, |r3| obtained in 2+1 
flavour dynamical staggered QCD on A^t=4 lattices. 

Figure ^ shows the phase diagram on the fi-T plane in physical units, thus Tc as a 
function of /i^, the baryonic chemical potential (which is three times larger then the quark 
chemical potential). The analysis is consisted of three steps. First one determines the 
transition points as a function of /x. Then the V ^ oo behaviour is inspected to separate 
the crossover and the flrst order phase transition region. Finally one transforms lattice 
units into physical ones. In physical units the endpoint [ |[ is at = 160 ±3.5 MeV, 
fXE = 725 ± 35 MeV. At /iB=0 we obtained = 172 ± 3 MeV. Note, that due to CPU- 
limitations in this analysis the light quark masses are approximately three times larger 
than their physical value and the lattice spacing is ^0.28 fm. Clearly more work is needed 
to extrapolate to the thermodynamic, chiral and continuum limits. 
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Figure 4. The T-fi diagram. Direct results are given with errorbars. Dotted line at small 
/i shows the crossover, solid line at larger fi the first order transition. The box gives the 
uncertainties of the endpoint. 
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Figure 5. Left panel: p/T^ as a function of T/Tc at yU = The continuum SB limit is 
also shown. Right panel: Ap = p{fi 7^ 0, T) — p{fi = 0, T) normalized by as a function 
of T/Tc for /iB=100,210,330,410 MeV and /iB=530 MeV (from bottom to top). 



The equation of state at T 7^ and yU 7^ is also determined [ lo\- In order to help 



the continuum interpretation of the figures we normalize the raw lattice results with the 
dominant correction factors between Nt=4: and the continuum in the T ^ 00 (Stefan- 



Boltzmann; SB) case (see Ref. [ IT^). Therefore, the results presented on the figures might 



be interpreted as continuum estimates and could be directly used in phenomenological 
applications. 

Figure ^ shows the pressure (p) at /i=0 and Ap{T ,fi)=p(T ,fi)-p(T ,^=0) normalized 
by T^. Note, that normalizing Ap{T,fi) by Ap{T 00,/i) leads to an almost universal 
/^-independent function [ 0. The left panel of Figure ^ shows e-3p normalized by T"^, 
which tends to zero for large T (e is the energy density). The right panel gives the baryon 
number density as a function of T/Tc for different fi-s. As it can be seen the densities 
exceed the nuclear density by up to an order of magnitude. 
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Figure 6. Left panel: (e — 3p)/T^ at /iB=0,210,410 MeV and /iB=530 MeV as a function 
of T/Tc (from bottom to top). Right panel: The baryon number density normalized by 
as a function of T/T^ for /i^ =100,210,330,410 MeV and /iij=530 MeV. As a reference 
value the line starting in the left upper corner indicates the nuclear density. 




6. Overlap improving multi-parameter reweighting: Taylor expansion 



This section summarizes some results of the Bielefeld- Swansea group [ |^, Q obtained 
in dynamical staggered QCD with p4 action on Nt=4 lattices. 

Instead of evaluating the determinants in eq. (1) explicitely one can approximate them 
by using a Taylor expansion as given by eq. (||). One calculates observables as a function 
of fi using eq. (||). The left panel of Figure ^ shows the chiral susceptibility as a function 
of the gauge coupling for three different /i values. Note, that the peak indicates the 
transition and that its position moves as /i changes. Since the gauge coupling is directly 
connected to T one can convert the position of the peaks into physical units and obtain 
the Tc as a function of fi (right panel of Figure |^). The authors of Ref. concludes that 
the curvature of the phase diagram at /x=0 is in good agreement with that of Ref. [ H 
(they show the endpoint of Ref. [ || by a diamond signaled as Fodor & Katz). 

The equation of state was also studied at small values. At the RHIC point both e 
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Figure 8. rricrit/T as a function of fi'^/T^. The band indicates the uncertainties. 

and p increase by around 1% from its fi=0 value. Along the transition line their change 
is consistent with zero within the current precision [ |^ . 

In a hypothetical QCD with degenerate quark masses and m7r~mx~190 MeV the 
critical endpoint is at fi=0. Note, that fi^O values correspond to larger critical masses 
{'merit)- A preliminary result is shown in Figure ||. Using the n dependence of merit one 
could, in principle, estimate the location of the endpoint for physical quark masses. 



7. Analytic continuation from imaginary /i 

This section summarizes the results obtained by the analytic continuation technique. 
Nt=4 lattices are used with 2 flavour [ |I^] and 4 flavour [|T^ dynamical staggered QCD. 



In Ref. [ |T0[ an alternative approach was developed, avoiding reweighting in alto- 
gether. This is achieved by simulating with imaginary fi, where there is no sign problem 
and hence no need for reweighting. In this case one may fit the nonperturbative data of 
an observable (even of the critical line) by truncated Taylor series in fi/T (see left panel 
of Figure pT). In the absence of nonanalyticities, the series may be analytically continued 



to real values of /i (right panel of Figure ^). For the transition temperature [ 10| one gets 
Tc(/i) = Tc(0) — 0.0056/i^/T (note, that similar value was obtained in four flavour QCD 
by Ref. [ |ll]). These results are in complete agreement with those of the previous two 
sections. 

Performing the analytic continuation for different volumes, a finite volume scaling anal- 
ysis should reveal the nature of the transition. In particular, it might then be possible to 
locate the critical endpoint by this technique, too. 



8. Conclusions 

Due to the notorious sign problem this is the first time that basic thermodynamic 
quantities could be obtained on the lattice at finite chemical potentials. Completely 
different new techniques lead to very similar result, which indicates that in the near 
future lattice QCD could be a major contributor to the field. Clearly, more work has to 
be done in order to reach the thermodynamic, chiral and continuum limits. 
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from the Bielefeld-Swansea group are also acknowledged [ |12| . 
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Figure 9. Left panel: Transition coupling for imaginary /i with different fits. Right 
panel: Transition temperature in physical units for fiy^O. 
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